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ABSTRACT 



We present the first detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL 
and RXTE. 4U 1909+07 is detected in the ISGRI 20-40 keV energy band with an average countrate of 2.6 cts sec~' . The pulse period 
of ~604 sec is not stable, but changing erratically on timescales of years. The pulse profile is strongly energy dependent: it shows 
a double peaked structure at low energies, the secondary pulse decreases rapidly with increasing energy and above 20 keV only the 
primary pulse is visible. This evolution is consistent between PCA, HEXTE, and ISGRI. The phase averaged spectrum can be well 
described by the sum of a photoabsorbed power law with a cutoff at high energies and a blackbody component. To investigate the pulse 
profile, we performed phase resolved spectral analysis. We find that the changing spectrum can be best described with a variation of 
the folding energy. We rule out a correlation between the black body component and the continuum variation and discuss possible 
accretion geometries. 

Key words, stars: neutron (4U 1909+07) - X-rays: binaries - Accretion 



1. Introduction 

The accreting neutron star X-ray binary 4U 1909+07 was 
first mentioned in th e 3 rd Uhuru catalog as 3U 1912+07 
dGiacconi et all 1 19741) . The position and name of the source 
were later refined t o 4U 1909+07 in the 4 th Uhuru catalog 
(iForman et al 1 119781) . Many other X -ray missions, such as OSO 
7, Ariel 5, and EXOSAT, among others, als o detected a source 
close to these coordinates. IWen et alj (l2000h showed that these 
detections are likely all originating from the same source and re- 
fer to it as X1908+075. We will use the name 4U 1909+07, to 
honor the original discovery. 

Despite the regular detections, the first dedicate d paper to 
discus s 4U 1909+07 was published in 2000 when IWen et all 
d2000l) analyzed RXTE-ASM data and found a stable period 
of 4.4 days. This period is interpreted as the orbital period of 
a binary system. Due to the high photoabsorption, however, 
no optical counterp art could be identified. In RXTE-PCA data, 
iLevine et all (|2004) found a second, shorter period of ~605 s in 
the X-ray flux, explained as the pulse period of a slowly rotating 
neutron star. Using Doppler delay curves, they could also refine 
the binary orbit parameters and estimated the mass of the com- 
panion star to be M* = 9-3 1 M and the radius of the compan- 
ion star to be < 22 R Q , adopting a canonical mass of 1 .4 Mp 
for the neutron star. One year later, iMorel & Grosdidierl <l2005h 
detected an OB star in the near infrared at the location of the 
X-ray source, thus confirming that the system is a High Mass X- 
ray Binary ( HMXB). The distance of the system was estimated 
to b e 7kpc dMorel& Grosdidierl l2005h . Prior to this discov- 
ery, ILevine et aL (2004) argued that the companion star could 
be a Wolf-Rayet star, which would make this system a possible 
progenitor system to a neutron star-black hole binary. With the 
identification of the companion, however, this intriguing possi- 
bility can be ruled out. The X-ray luminosity L of 4U 1909+07 is 



around 2.8 x 10 J6 ergs _1 for 4.5-200 keV. Although the system 
shows n o eclipse, the X-ray flux is still strongly or bital phase de- 
pendent dLevine et all I2004T) . ILevine et all d2004l) also analyzed 
orbital phase resolved PCA spectta and found that around orbital 
phase tpofb = 1 the photoabsorption increases by a factor of 2 or 
more to A^h > 30 x 10 22 cm~ 2 , explaining the decreased ASM 
flux. This increase in absorption can be very well described by 
a spherical wind model and an inclination of 54° < i < 70°, 
depending on the parameters of the wind model. 

In archival Chandra data, Torrejon et al] d2010l) found that 
the soft energy spectrum of 4U 1909+07 shows evidence for a 
Compton shoulder on the soft energy edge of the iron Ka line. 
This feature has so far only been seen in GX 301 -2 but in no 
other HMXB with Chandra. iTorreion et aT] <l2010h showed that 
the Compton shoulder is consistent with the X-ray source being 
embedded in a Compton thick medium, which is also responsi- 
ble for the observed photo absorption. Additionally these authors 
have shown that 4U 1909+07 is an exception from the correla- 
tion between A^h and t he equivalen t width of the iron line found 
in most other HXMB dlnouel 1985). The origin of this excep- 
tion is still a mystery and can only be investigated with more 
high-resolution spectra. 

Although there was renewed interest in 4U 1909+07 in re- 
cent years, many basic parameters are still unknown. For ex- 
ample, despite the detection of the pulse period, its evolution 
with time was not yet studied, even thoug h it can give insigh t 
on the prevailing accretion mechanism (Bildst en et all |_1997). 
Additionally, no detailed analysis of the high energy spectrum 
has been carried out so far. RXTE data provide enough statis- 
tics to even perform pulse phase resolved spectroscopy, allowing 
us to obtain a better understanding of the accretion region and 
mechanism together with the energy dependence of the pulse 
profile. In this article we are aiming at improving our knowl- 
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Table 1. Observation log of the RXTE and INTEGRAL data 
used. The dates given are the start dates of the respective obser- 
vation. 



Obs-ID / ScW 




Start date 


Exposure 


RXTE 








70083-01-01-00 


-70083-01-28-00 


2002-12-23 


105ksec 


70083-02-01-00 


-70083-02-19-00 


2003-01-30 


75 ksec 


INTEGRAL 








004800030010 - 


007000670010 


2003-03-06 


1352 ksec 


017200250010 - 


020200020010 


2004-03-11 


1192 ksec 


022600370010 - 


025900180010 


2004-08-20 


966 ksec 


029100170010- 


032300140010 


2005-03-02 


340 ksec 


036100020010- 


038300520010 


2005-09-27 


880 ksec 


047600030010 - 


050600200010 


2006-09-06 


397 ksec 


053700090010 - 


055600710010 


2007-03-07 


608 ksec 


060100090010 - 


062300080010 


2007-09-15 


967 ksec 


072200090010 - 


074300080010 


2008-09-11 


1072 ksec 


078200500010 - 


080400080010 


2009-03-10 


690 ksec 



edge in these points, using data from INTEGRAL and RXTE. 
In Sect. |2]we present the data and reduction methods. The pulse 
period evolution and the pulse profiles at different energies are 
analyzed in Sect. [3] In Sect. [4] we perform phase averaged and 
phase resolved spectroscopy. We summarize and discuss our re- 
sults in Sect. [5] 



2. Observations and data reduction 

For our study of the 4U 1909+07 system we used all available 
public data from the X-ray missions INTERnat ional Gamma 
Ray 



Astrophysics Laboratory (INTEGRAL, Winkler et al 
2003) and Rossi X-Ray Timing Explorer (RXTE, [Bradt et al 



1993). One of the main detectors of INTEGRAL is ISGRI, a 
coded mask instrument sensitive in the 15 keV - 1 MeV energy 
range and part of th e Imager on Board t he Int egral Satellite 
(IBIS; lUbertini et all l2003b iLebrun et all 120031) . IBIS/ISGRI 
was the first instrument to produce high resolution images 
of the X-ray sky above 20keV. Even though no pointed ob- 
servations on 4U 1909+07 were performed with INTEGRAL, 
thanks to IBIS/ISGRI's large field of view of almost 30° X 30°, 
there are more than 8.5 Msec of off-axis data available. A 
huge part of these data were obtained in the core program on 
the Galactic Center performed in the early Announcements of 
Opportunity (AOs; Winkler, 2001b . Data from later AOs were 
mainly taken from the exten ded monitoring campaig n of the mi- 
croquasars GRS 1915+10 5 ([Rodriguez et all 12008b and SS 433 
dKrivoshevev et all 120091) . See Table Q] for a detailed list of the 
data used. 

An IBIS/ISGRI mosaic showing the hard X-ray sky in spring 
2003 for the region around 4U 1909+07 is shown in Fig. [TJ 
Twelve sources with a detection significance > 6cr were found in 
that mosaic, with 4U 1909+07 being the fourth brightest source 
with ~2.6cts sec -1 . The data were extracted using the standard 
pipeline of the Off-line Scientific Analysis software (OSA) 7.0 
for spectra and images and using ii .light to obtain lightcurves 
with a higher temporal resolution. We extracted lightcurves with 
a resolution of 20 sec as a trade-off between good signal-to- 
noise ratio and good time resolution to measure the pulse pe- 
riod. We shifted the time to the barycenter of the solar system 
and corrected for th e orbital motion of th e neutron star, using the 
ephemeris given bv lLevine et all (|2004|) . 

For phase resolved spectroscopy we used RXTE observa- 
tions carried out in early 2003 which were previously ana- 



-1.1 -0.82 -0.53 -0.24 0.051 0.34 0.63 0.92 1.2 

Fig. 1. ISGRI image of the X-ray sky in the 20^-0 keV energy 
band, averaged over 353 ScWs from Revs. 51-68. Coordinates 
are in galactic longitude and latitude. To avoid confusion only 
the 6 brightest of overall 12 detected sources are labeled. In 
this image 4U 1909+07 is the fourth brightest source in the field 
of view after GRS 1915+105, the transient 4U 1901+03, and 
H 1907+097. 



lyzed bv iLevine et a D d2004). These data had a total exposure 
of 196 ksec. The only other RXTE data are from a shorter ob- 
servation in 2000, which are not suitable for our study. Both 
major instruments of RXTE were fully operational at the time 
of the observations in 2 003, the Proportional Counter Array 
(PCA; iJahoda et ail 120061) . sensitive in the 2-60 keV energy 
band and the High Ene rgy X-Ray Timing Experiment (HEXTE; 
iRothschild et all 1 1998th sensitive between 1 5-250 keV. Both in- 
struments have a large effective area and a high time resolu- 
tion, but no imaging capabilities. With a field of view of ~1° 
no other sources are visible to RXTE's instrument while point- 
ing at 4U 1909+07. RXTE data were reduced using the standard 
HEASARC software (v. 6.6), obtaining lightcurves of the PCA 
with 1 sec resolution and pulse phase resolved spectra for PCA 
and HEXTE in the 4.5-150keV energy range. The lightcurves 
were shifted to the barycenter of the solar system and corrected 
for the neutron star's orbital motion in the same way as the 
INTEGRAL data. 



3. Timing Analysis 

3. 1 . Pulse period evolution 

Monitoring the evolution of the pulse period allows us to deter- 
mine the type of accretion, i.e., if direct wind accretion prevails 
or if a stable accretion disk is formed. For 4U 1909+07 the his- 
tory of pulse period me asurements starts on ly in 2000 with its 
detection in RXTE data (ILevine etali r2004) and is continued 2 
years later by another measurement o f RXTE. At both ti mes a 
period of P =s 604.7 s was measured (ILevine et al ., 2004- The 
archival INTEGRAL data provide an optimal basis to perform 
many more measurements of the pulse period between 2003 and 
2008 and to follow its evolution more closely. To do so, we 
split the IBIS/ISGRI lightcurve into segments between 300 ksec 
and 800 ksec length, providing a good balance between accurate 
pulse period determinati on and tempor al resolution, and used the 
epoch folding technique (iLeahvl, 0*987) with 24 phase bins to de- 
termine the pulse period for each of these segments individually. 



F. Furst et al.: 4U 1909+07: a well-hidden pearl 



3 




Time in sec since MJD 52704 

Fig. 2. Small part of the lightcurve in the 2CM-0 keV energy band 
as measured with IBIS/ISGRI. Data are shown as black stars 
with gray bars indicating the uncertainties. Time resolution is 
20 sec. The line shows the folded pulse profile over the whole 
data set. 
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Fig. 3. Evolution of the pulse peri od from 2001 until early 2009. 
The historic RXTE data points of lLevine et alJ ([2004) are shown 
as red crosses, the new measurements obtained with INTEGRAL 
are shown as stars. 



The counting statistics in ISGR1 smear out the single pulses 
in the lightcurve and make them invisible to the naked eye, see, 
e.g., Fig. |2 Only after folding «500-1000 pulse periods, a stable 
pulse profile emerges which allows for a reliable determination 
of the pulse period. This profile is remarkably stable over our 
data set and not time dependent. 

The pulse period evolution between 2001 and 2007 is shown 
in Fig. |3] The uncertainties of the pulse period were estimated 
by fitting a Gaussian distribution to the x 2 distribution of the 
epoch folding results. The width of the best-fit Gaussian is a 
good estimate for the uncertainty of the measurement. 

For the pulse period determination we assumed the pulse 
period to be constant during each individual segment. This as- 
sumption is justified, as the average P in our data between two 
consecutive data points is on ly f « 1.3 X 10~ 8 sec sec 1 , very 
close to the value given by iLevine et alJ (2004). With the in- 
dividual lightcurves being on average shorter than 50 days, the 
change between start and end of a lightcurve would be only 
AP x 0.05 sec. This change is on the order of the uncertainty we 
obtained from the Gaussian fits. Even when assuming a larger 
value of P = 1 sec yr~' =3.17 sec sec -1 the obtained pulse pro- 
files did not change significantly in shape, making a distinction 
impossible. 

The cros s es in F ig. [3] are RXTE measurements taken from 
ILevine efa l. (2004) and indicate an almost unchanged pulse 
period between January 2001 and December 2002. Starting 
with the INTEGRAL data, shown as stars in Fig. [3] a spin- 
up trend between March 2003 and September 2005 is visible 
from P « 604.7 sec to to P » 603.9 sec. This trend is bro- 
ken after September 2005 and the pulse period varies around 
P = 604.0 sec. These changing trends make the overall behavior 
consistent with a random walk . To confirm this observ ation, we 
used the algorithm proposed bv lde Kool & Anzen(ll993l) . It eval- 
uates the relative change in pulse period over different time in- 
tervals St between single measurements. To simplify the calcula- 
tions the pulse period is converted to angular velocity u> = 2n/P 
and its change is expressed as the average value of absolute dif- 
ferences: 

5(o(5t) = (\a>(t + St) - oj(t)\) 

Each pair (a>(t + St),io{t)) is additionally weighted in or- 
der to account for the data being grouped in clusters, see 



Ide Kool & Anzerl (1 19931) for a detailed description. The results 
of this analysis are shown in Fig. [4] 

If the data are compatible with a random walk, they will fol- 
low a straight line with a slope of 0.5 in log Sa> - log St space. To 
test if the data are consistent with such a straight line, we fitted 
a function of the form 

\og{8a){8t)) = A + 0.5 • log(<5f) 

where the free parameter A describes the noise level of the data. 
The best fit, superimposed in Fig. [4] has a Pearson's correlation 
coefficient of r ~ 0.95 and clearly shows that the pulse period 
evolution is consistent with a random walk. It gives a noise level 
of A = —9.1, fully in agreement with the expected value for a 
source with a luminosity of L « 2.8xl0 36 ergs " 1 , assuming wind 
accre tion with a turbulent accretion wake dde Kool & Anzerl 
1 19931 and references therein). 

The uncertainties of the £w-values are originating not only 
in the determination of the pulse period but also in the uneven 
and coarse sampling of the pulse period evolution. To estimate 
the latter effect, we chose a Monte Carlo approach and simu- 
lated many mock-up pulse period evolutions which followed a 
perfect random walk and sampled these with the same rate as 
the real data (see |de Kool & Anzerl fl993). The standard devia- 
tion in each St bin from these simulations gives an estimate of 
the uncertainty in each St bin of the data and is roughly of the 
same order of magnitude as the uncertainties from the pulse pe- 
riod determination. 

3.2. Pulse profiles 

ILevine et a I1 (l2004l) showed in their Fig. 2 that the average pulse 
profile between 3.7 - 17keV shows two distinct peaks, with the 
second one being slightly broader than the first one and exhibit- 
ing a complex shape with two subpeaks. With higher energy res- 
olution a remarkable and strong energy dependence becomes ev- 
ident, clearly seen in pulse profiles extracted from RXTE/PCA 
data with a sampling period of P = 604.685 s and 32 phase bins. 
Thanks to the large effective area of the instrument, very good 
pulse profiles with a high signal-to-noise ratio could be extracted 
in 30 narrow energy bands. Three examples are shown in the up- 
per panels of Fig. [5] At very low energies below 5 keV a broad 
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Fig. 4. Plot of t he pulse period evo lution in log Sa> - log 5t space, 
as proposed bvlde Kool&Anzen(fl993l An ideal random walk 
process would show a slope of 0.5 in this plot. The data follow 
this slope closely, as shown by the red straight line with a slope 
of 0.5. 

plateau is visible, with increasing intensity at later pulse phases 
(Fig. |5^). With increasing energy this plateau separates clearly 
into two peaks, while the second peak is getting dimmer com- 
pared to the first one (Fig. [5J)). The relative power of the sec- 
ondary peak keeps declining up to ~20keV (Fig. [5};) and it is 
not visible anymore above 20keV (Fig.[5ji). Despite the strong 
evolution with energy of the pulse profile, the deep minimum 
around phase <p = 0.3 is not energy dependent. The IBIS/ISGRI 
pulse profile, which was extracted from the average data of the 
first 100 days of measurement in 2003, is shown in Fig.|5}l for 
comparison. We used the same epoch as for the PCA analysis, 
but a period of P — 604.747 s, as determined in our analysis of 
the data. The pulse profile at energies above 20 keV is consistent 
between IBIS/ISGRI and RXTE/HEXTE. 

To make the shape transition of the pulse profile at low and 
medium energies better visible, we plotted them in a color coded 
map (see Fig. |6j. The x-axis shows the pulse phase, the log- 
arithmically scaled y-axis displays the energy in keV and the 
mean count rate is color-coded. The extensions of the pixels in 
y-direction represent the energy band of the given pulse profile. 
A smooth transition from the broad, two peaked structure to the 
single peaked structure is clearly seen. The minimum around 
phase <p = 0.3 in black stays very stable. 



4. Spectral analysis 

4. 1 . Phase averaged spectrum 

The pulse phase averaged spectrum of 4U 1909+07 shows a 
shape similar to many other neutron star sources, consisting 
of a powerlaw continuum attenuated by photoabsorption at 
low energies and a turnover at high energies. For this kind 
of continuum various phenomenological models exist, all de- 
scribing a similar shape. Among the most widely used mod- 
els are the cutoff-powerlaw cutof fpl, with a smooth turnover 
at higher energies (determined by the folding energy £f id), 
the highecut mod el, where the tur nover begins suddenly at a 
given energy E cut (IWhite et al.LI 1 983l a nd references there i n) and 
the Comptonization model compTT (IHua & TitarchuM [1995) 
with more free parameters to describe the tem perature and 
optical thickness of the Comptonizing plasma. iLevine et al.l 




Phase 

Fig. 5. Energy resolved pulse profiles with RXTE PCA (a-c) 
and RXTE HEXTE (d). The profiles are shown twice for clar- 
ity. In panel d) the INTEGRAL pulse profile in the 20-40 keV 
energy range is shown in gray for comparison. Note that the 
INTEGRAL profile is shifted by hand to match the peak in the 
RXTE profile, but that the profiles are not phase aligned. 



(2004) modeled the spectra using a cutoff powerlaw and a 
bremsstrahlung model, which incorporates only the temperature 
to describe the spectral shape. They found that both models fit 
the data equally well. We applied all four models, cutof fpl, 
highecut, compTT, and bremss to the RXTE and INTEGRAL 
data simultaneously, using ISIS 1.6.0-3 and discarding all ob- 
servations between orbital ph ase 0.88 < foyb < .12, where the 
TYh is dramatically increased (ILevine et all l2004). The photoab- 
sorption at low energies was modeled u sing a revised versi on of 
the tbabs mode Q with abundances b y Iwilms et all (l2000h and 
cross-sections bv IVerner et al.1 dl996h . A strong iron line close 
to the ene rgy of neutral iron o f 6.4 keV is also presented in the 
spectrum (ILevine et all [2004), and was modeled for all contin- 
uum models using a Gaussian emission line. 

All continuum models described the data very well, except 
the bremsstrahlung model, which only gave a best fit reduced 
X 2 value of x 2 ed = 1-75. We ascribe this to the fact that the 
hard X-ray spectrum of 4U 1909+07 can not be modeled accu- 
rately with the bremss model, noticeable only thanks to the im- 
proved statistics with the additional IBIS/I SGRI data. These data 
were not available to iLevine et a D d2004). The best fit parame- 
ters and^ 2 d values of the other models are presented in Table [2] 
and Fig. [7J1 shows the RXTE and IBIS/ISGRI spectra together 
with the best fit cutof fpl. 

The x 2 ed values are below 1.5 for all models presented in 
Table [2] indicating a statistically acceptable fit. A look at the 
residuals in Fig. |7j), however, shows that the models do not de- 
scribe the data satisfactorily below lOkeV. These residuals can 
be eliminated by adding a blackbody component (Fig. [TJ; and 



1 see http: //pulsar. sternwarte.uni-erlangen.de/wilms/research 
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Pulse phase 



Fig. 6. Color coded pulse profile map of PCA and HEXTE data. 
The color indicates the mean countrate in each bin, ranging from 
black (lowest) to white (highest). The two highest energy bins 
are from HEXTE data. The four horizontal blue lines indicate 
the energies for which the pulse profiles are shown in Fig. [5] 
The map shows one pulse only. 

Table. [2]). The cutoffpl and the highecut model give similar 
blackbody temperatures between kT «1.4~1.7keV, whereas in 
the compTT the temperature increased to kT — 3.1 + 0.4 keV. As 
the ;^ 2 of the compTT + bbody model is not as good as the x 2 of 
the other two models, a temperature around kT as 1 .5 keV seems 
most realistic. 

4.2. Phase resolved spectra 

Due to the rotation of the neutron star, the accretion columns and 
hot-spots are seen under a constantly changing viewing angle. 
Since the physical conditions can not be assumed to be constant 
throughout the emission regions, the X-ray spectrum can change 
significantly with pulse phase. These dramatic changes are indi- 
cated in the energy dependence of the pulse profile (Fig.[5]and[6jf. 
Pulse phase resolved spectroscopy allows us to disentangle the 
contributions from the di fferent physical components (see, e.g., 
iRrevkenbohm et all, |2004|) . To perform phase resolved spectral 
analysis, we divided the RXTE/PCA and RXTE/HEXTE data 
into 7 phase bins, as indicated in Fig. [8^. The IBIS/ISGRI data 
did not provide good enough statistics for high resolution phase 
resolved spectroscopy. A cutoffpl, a highecut, and a compTT 
model, each improved by photoabsorption, an iron line and a 
black body were fitted to all phase bins. All three models showed 
similar behavior as function of pulse phase. In terms of x 2 ed > 
however, the highecut failed in most phase bins very clearly, 
even though it gave the best ^ 2 ed values for the phase averaged 
spectrum. A possible explanation for this effect is that the phase 
averaged spectrum is a superposition of the different spectra at 
different phases and that this superposition is best described with 
the highecut model. To check this assumption, we simulated a 
spectrum consisting of a superposition of the best fit cutoffpl 
models for the individual phasebins, each parameter weighted 
by the respective count rate in the respective phase bin. The 
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Fig. 7. a) Combined RXTE PCA (red) and HEXTE (green) and 
INTEGRAL ISGRI (blue) spectra. The best fit cutof fpl-model 
with a blackbody component is shown in black, b) Residuals of 
the best fit model without a black body, c) Residuals of the best 
fit model with a black body. 



emerging spectrum could be described significantly better with 
a highecut model 0f 2 ed = 0.72) than with the cutoffpl model 
(y 2 ed = 0.82). This result shows that the superposition of the sim- 
pler cutoffpl models introduces new features in the phase av- 
eraged spectrum, which are not source intrinsic. The cutoffpl 
model gave the best fits over all phasebins and only those results 
are presented here. 

Describing the spectra without a blackbody component re- 
sulted in good x 2 ed values only for phase bins A, B, C, and G, but 
returned formally unacceptable x 2 ed values for phase bins D, E, 
and F, i.e., for spectra during the secondary peak. As discussed 
in the previous section the phase averaged spectrum can also be 
fitted relatively well without a black body component, but the in- 
clusion of such a component improves the y 2 value. This effect 
can now be explained by viewing the phase averaged spectrum 
as a superposition of spectra with and without a prominent black 
body component, and consequently showing only comparatively 
weak evidence for this component. Phase resolved analysis dis- 
entangles this superposition and gives more information on the 
real physical behavior of the source. 

A close analysis of all spectra revealed that the folding 
energy and the photon index of the cutoffpl model depend 
strongly on each other and that the statistics of the spectra do 
not allow us to separate the influence of these two parameters. 
We consequently carried out two independent fits, one with the 
folding energy frozen to its value in the phase averaged fit and 
one with a frozen photon index, also set to the value of the phase 
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Table 2. Fit parameters for the phase averaged spectra with different models. No significant width of the iron line could be measured 
for the highecut + bbody and the compTT model and Fe <x was consequently frozen to 0.1. Not shown are the cross calibration 
factors and the renormalization of the PCA background, which are similar for all models and close to 1 . 



Model 


cutofFpl 


cutofFpl 


highecut 


highecut 


compTT 


compTT 


parameter 

a t rinll -2-1 




+bbody 




+bbody 


in T + l.K 


+bbody 



VH I 



1 .9 



Accm [ph@lkeV] 0.084^ 0.029^ 0.057+°™ 0.036^ (0.54 ± 0.03) x 10~ 2 (4.8 ± 0.4) x 10~ : 

T 1.59 ±0.07 1.15 ±0.12 1.565* l36+ -o\l 

E foM [keV] 28 ± 3 20 ± 2 27 ± 3 23 ± 3 

£cu,off[keV] - - 7.8 ±0.5 7.5 + ^ 

r [keV] - - 1.19 ±0.04 1.15 ±0.06 

kT[ke.V] - - 9.7+";; 10.1 ±0.3 

r - 3 01 + ° 13 3 00 +o u 

Fecr[keV] 0.9^ 0.19!°)^ 0.41 ±0.17 0.1 0.1 0.15™^ 

Fe Energy [keV] 6.76 ±6". 15 6.47 ± 0.07 6.43^ 6.42+° '° 6.44+° ^ 6.44 ± 0.06 

Fe A [ph/s/cm 2 ] (0.74+" Jf ) x 10~ 3 (4.3^ 2 ) x 10~ 4 (0.6+g j) x 10" 3 (3.9+^) x lO" 4 (3.4^) x 10" 4 (4.5+^) x 10" 4 

bbody kT [keV] - 1.65+°^ - - 3.1 ±0.4 

bbody norm - (0.71+g-j^) x 10 3 - (5 ± 2) x lO" 4 - (3.7+| 4 ) x lO" 4 

A^/dofOtL) 147.4/100 (1.47) 94.6/98 (0.97) 102.7/99 (1.04) 85.4/98 (0.87) 124.5/100(1.25) 98.8/97 (1.02) 



averaged fit. Both models gave very good results in terms of 
X^,. To further reduce the number of variable parameters, we 
also froze the energy of the Gaussian iron line to 6.47 keV, the 
value of the phase averaged fit. The iron line energy was con- 
sistent with that value in all phase bins when allowed to vary. 
Overall the free parameters of the fit were: the normalization of 
the cutoffpl, the column of the photo absorption, the normal- 
ization and temperature of the black body, the normalization and 
width of the iron line and either the folding energy or the photon 
index of the cutoffpl. 

The evolution of the spectral parameters is shown in Fig. [8] 
for a frozen folding energy and in Fig. [9] for a frozen photon 
index. It is clearly seen that for both models neither the photoab- 
sorption Nu nor the power of the iron line seem to vary much 
over the pulse phase. The strong spectral changes indicated in the 
energy resolved pulse profiles are thus mainly due to a change in 
the continuum of the cutoffpl model, either due to an increase 
in the photon index (Fig. [8}i) or due to a shift of the folding 
energy to higher energies (Fig. [9]d). The folding energy is high- 
est in the primary peak during phase bins A and B, making its 
spectrum distinctly harder than the rest of the pulse phase. In the 
secondary pulse the folding energy moves to values as low as 
~17keV, forcing a strong attenuation of the spectrum above this 
energy. 

The behavior of the black body component is also very inter- 
esting in both models. Figures [8] and [9] show the normalization 
of the black body component and its temperature (panels e and f, 
respectively). While the temperature varies insignificantly with 
pulse phase, the normalization, expressing the contribution of 
the black body component to the overall flux, shows a strong de- 
pendency on pulse phase with variation up to a factor of 3. The 
black body is strongest in phase bin D, i.e., during the beginning 
of the secondary pulse, and lowest in phase bin A, i.e., during the 
rise and maximum of the primary peak. As there is an indication 
for a broadening of the iron line during the phase bins where the 
black body is weakest, one might speculate about a compensa- 
tion of the blackbody component by the Gaussian. Confidence 
contours between the two parameter, however, showed no corre- 
lation. The indicated broadening of the iron line is more likely 
due to a change in the iron edge or an iron KJ3 line, which can not 
be resolved with RXTEfPCA. Using spectra with higher energy 
resolution a more in depth study of the iron line behavior with 
pulse phase will be possible. 



Regarding the normalization of the powerlaw, it clearly re- 
flects the pulse profile in the model where the photon index was 
frozen (see Fig. [9};). In the model with the frozen folding en- 
ergy, on the other hand, the correlation is not so clear. It rather 
seems that in this case a combination of the photon index and the 
power law normalization models not only the spectral change, 
but also the brightness variation. Because of this entanglement, 
the model with the variable folding energy seems the more real- 
istic description of the spectrum. 

To exclude strong dependencies we calculated confidence 
contours for the phase averaged spectrum for the blackbody nor- 
malization versus the powerlaw normalization as well as for the 
blackbody normalization versus the folding energy, based on the 
models with frozen folding energy and frozen powerlaw normal- 
ization, respectively. Figure[10]shows the confidence contour for 
the model with the variable photon index, i.e., the frozen fold- 
ing energy, together with the best fit value as a bright cross. The 
contours are indicating only a very weak correlation between 
the black body normalization and the photon index. In the back- 
ground the best fit value of the normalization of the Gaussian 
iron line is shown in gray-scale for every relevant point in the 
black body normalization - photon index space. The scale on 
the right of the plots gives its corresponding values. With in- 
creasing blackbody power and spectral softening a weaker iron 
line is necessary. Within the best fit contours, however, the iron 
line is not changing dramatically, thus ruling out a substitution of 
the blackbody with a stronger iron line for formally acceptable 
Rvalues. 

Additionally to the phase averaged contour, Fig. [10] shows 
the best fit values for all 7 phase bins as dark crosses, labeled 
with the respective letter. As can be clearly seen, in phase bins A 
and B the spectrum is hardest, while the black body component 
is weakest. Phase bins C-F are clustered together in a regime 
with a slightly softer spectrum and stronger black body than the 
phase averaged spectrum. Phase bin G, describing the minimum 
between two consecutive pulses shows only very weak portions 
of the black body and an intermediate spectral hardness. This 
distribution makes very clear that both, the photon index and the 
black body component change significantly over the pulse phase 
in a way which seems to indicate a correlation between them. 
This correlation is most likely of physical origin and not intro- 
duced due to improper modeling of the spectrum, as the confi- 
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Fig. 8. Parameters of the phase-resolved spectra fitted with the 
cutoffpl model and the folding energy frozen to Ef \d = 
20.0keV . a) Pulse profile in the 13.2-14.5keV energy range. 
The different shaded areas indicate the phase bins used to extract 
the spectra, b) Photo electric absorption column A^h, c) power 
law normalization, d) power law index Y,e) black body normal- 
ization, f) black body temperature, g) iron line normalization, 
h ) iron line width, and i) x 1 
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dence contour of the phase averaged spectrum showed no corre- 
lation between these parameters. 

Figure [TT1 shows the same plot as Fig.[10]but for the model 
with a variable folding energy, i.e., a frozen photon index. As for 
the variable photon index, the iron line does not change signifi- 
cantly within the acceptable^ 2 contours and is only a weak func- 
tion of both parameters. The best fit phase resolved values show 
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Fig. 9. Same as Fig. [8] but for the cutoffpl model with the 
photon index frozen toF = 1.15. a) Pulse profile in the 13.2- 
14.5 keV energy range, b) photo electric absorption column A^h, 
c) power law normalization, d) folding energy, e) black body 
normalization, /) black body temperature g) iron line normaliza- 
tion, h) iron line width, and iJ^jL-value. 



a similar distribution as in Fig. [10] It is clearly seen that phase 
bin A shows the hardest spectrum with a negligible black body 
component, while phase bin B moves a bit to a softer spectrum 
with a stronger black body. Phase bins C-F are again clustered 
together and phase bin G shows intermediate spectral hardness 
with a very weak black body. This distribution of the phase re- 
solved data indicates an anti-correlation between the black body 
normalization and the folding energy. As in the other model, this 
correlation is most likely intrinsic to the source as the phase av- 
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Fig. 10. Confidence contour and best fit value for the phase av- 
eraged spectrum between the photon index of the powerlaw and 
the normalization of the blackbody component, when fitted with 
a frozen folding energy. The contour lines show from the cen- 
ter outwards 68% (lcr), 90% and 99% confidence level. Gray- 
scaled is the normalization of the Gaussian iron line for the best- 
fit model at every given photon index - blackbody normaliza- 
tion point. The points labeled A through G are the best fit values 
of the respective phase bins and their uncertainties in the same 
space and for the same model. 




Fig. 11. Same as Fig.[10]but for the folding energy of the power 
law versus the normalization of the black body, while the photo 
index of the cutof fpl model was frozen. 

eraged confidence contours rule out a model dependent correla- 
tion. 



5. Discussion & Conclusion 

We have presented the first detailed study with INTEGRAL and 
RXTE of 4U 1909+07 and have shown that the pulse period evo- 
lution shows strong indication for a random walk like behavior. 
Such a behavior has been seen in other HMXB like Vela X-l and 
is a strong indicator f or a wind accreting so urce without a per- 
sistent accretion disk dGhosh & Lamblll979h . An accretion disk 
would provide a more constant transfer of angular momentum 
and thus a long-term spin-up or spin-down trend, as seen in other 



sources like 4U 1907+09 (iFritz et all 120061) . Accretion disks 
form to fulfill the conservation of angular momentum of the ac- 
creted wind, but in a system with a strongly magnetized neutron 
star the ionized matter will couple to the magnetic field lines be- 
fore a stable disk can form. This requires field strengths on the 
order of 10 I2 G, not unusual for this kind of system. In these 
strong magnetic fields cyclotron resonant scattering features 
(C RSFs) can form, seen in m any other sour ces, e.g., Vela X- 
1 dKrevkenbohm et all 119991), 4U 0115+63 dSantangelo etaO, 
119991) . or MXB 0656-072 dMcBride et alll2006l) . We found no 
evi dence for a CRSF in th e spectra of 4U 1909+07. As shown 
by ISchonherr e t al. (2007) this does not rule out a strong mag- 
netic field, as CRSFs can be filled by photon spawning and their 
shapes and depths depend strongly on the geometry of the ac- 
cretion column and on the viewing angle. Further investigations 
with high-resolution data will allow a more thorough search for 
a CRSF in 4U 1909+07. 

A strong magnetic field in 4U 1909+07 is supported by the 
strongly pulsed flux, as the X-ray flux is believed to be produced 
in a very small region above the hot-spots on the neutron star 
surface, where it is penetrated by the mag netic field lines an d 
the matter is accreted onto the neutron star (Lamb et al., 1973). 
In the accretion column above the hot-spots, the density and tem- 
peratures are highest and the most X-rays are produced through 
Comptonization of the therm al photons from the therm al mound 
on the neutron star surface ( Bec ker & Wolffl. l2Q05alTbl) . If the 
magnetic axis is not aligned with the rotational axis of the neu- 
tron stars, the hot spots and the accretion columns move in and 
out of sight of the observer, resulting in a pulsed X-ray flux. 

The exact geometry of the accretion column depends 
strongly on the way the matter couples to the magnetic 
field and can not be constrained from our measurements 
(Basko & SunvaevL 119761: iMeszarosl 11984 . Besides the sim- 
ple filled funnel other possible configurations include a hol- 
low or a partly hollow funnel. With an average luminosity of 
~2.8 x 10 36 ergs~' between 4.5-200keV the system can host 
accretion rates that are large enough for a shock to form in the 
accretion column, due to the radiation pressure from the material 
close to the neutron star surface. This shock prevents escape of 
the Comptonized photons parallel to the magnetic field lines, so 
that they can only emerge from the sides of the accretion column 
in a so-called "fan beam" (Meszaros, 1984). If the accretion rate 
is lower, no shock will form and the radiation can escape along 
the accretion column in the "pencil beam". 

5. 1 . Classical picture 

In a classical assumption the two distinct peaks in the pulse 
profile at low energies hint at the fact that accretion happens 
onto both magnetic poles of the neutron star, but under different 
physical conditions. These differences can explain the different 
spectra in the two peaks. The secondary peak shows a distinctly 
lower folding energy which corresponds roughly to the tempera- 
ture of the electron gas in the accretion column. A misalignment 
between the magnetic center and the gravitational center of the 
neutron star could lead to different sizes of the accretion column 
and thus to different densities and temperatures. 

We have also seen that the blackbody component of the spec- 
trum of 4U 1909+07 is phase shifted with respect to the two 
peaks in the pulse profiles. It is likely that the blackbody origi- 
nates from the thermal photons of the thermal mound at the hot 
spot. Through the rotation of the neutron star, the hot spot is 
only visible at certain pulse phases. Assuming a hollow funnel 
accretion geometry, with a large enough accretion rate to form a 
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shock, we can argue that we see the thermal mound only at the 
phases were we look along the hollow accretion column, while 
the most hard X-ray flux is only visible when we see the sides 
of the accretion column, due to the "fan beam" geometry. As we 
averaged every phase bin over a large range of time, we can nei- 
ther rule out nor investigate a change of the accretion geometry. 

5.2. Strongly relativistic picture 

Caution has to be taken, however, with these simple classical 
models as the behavior of plasma in magnetic fields on the or- 
der of B « 10 12 G is still not well understood and all models 
are based on strong simplifications. Additionally gravitational 
light bending must be ta ken into account when analyzing pulse 
profil es and their origin dKraus et all 12003b iMeszaro s & Nagel, 
Il98l . Close to the surface of a canonical neutron star of mass 
M = 1.4 Mo and radius R = 10km, i.e., a radius only 2.4 
times the Schwarzschild radius, light bending leads to a visi- 
ble surface of the neutron star of 83%. This huge amount of al- 
ways visible surface reduces the pulsed flux from the two hot- 
spots. Assuming a "fan-beam" emission characteristic of the 
accretion column, gravitational light bending can increase the 
pulsed flux dramatically, especially at phases where one col- 
umn would be hidden behind the neutron star in the classical 
picture dRiffert & Meszarosl Il988h . In a very special configu- 
ration where the rotational axis takes an angle of 45° to both 
the observer and the magnetic axis, pulse profiles very simi- 
lar to the high energetic pulse profile of 4U 1909+07 (Fig. |5Ji) 
can emerge, in which only one sharp p eak is visible (com- 
pare Fig. 5 of iRiffert & Meszarosl Il988h . In the simple mod- 
els of Riffert & Meszaros ( 1988) the emission characteristic of 
the accretion column is taken to be isotropic and indepen- 
dent of energy, so that a direct translation to 4U 1909+07 is 
not easily made. Attempts to decompose the puls e profile for 
other s ources were made bv lKraus et alJ (120031) and lSasaki et alJ 
(l20Toh . among others. These decompositions prove difficult due 
to the many free parameters, like emission pattern, geometry 
of the accretion column, configuration of the magnetic field, 
and equation of state which can not be closely constrained. 
Nonetheless the sharp peaked pulse profile of 4U 1909+07 above 
20keV is a hint that the hard X-ray emission is produced in a 
"fan-beam" accretion column and that during the rotation of the 
neutron star one column vanishes behind the star and its emis- 
sion is bent around and focused due to the strong gravitation. 

More detailed models for pulse decomposition will be ap- 
plied to 4U 1909+07 in the future, together with data from other 
X-ray satellites that provide a higher spectral resolution, such 
as Suzaku. With these data a detailed investigation of the iron 
line complex and the Compton shoulder will also be possible 
and lead to more insight on the physical state of the absorbing 
medium. 
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